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ABSTRACT 

Context. HD 150136 is a triple hierarchical system and a non-thermal radio emitter. It is formed by an 03-3.5 V + 05.5-6 V close 
binary and a more distant 06.5-7 V tertiary. So far, only the inner orbital properties have been reliably constrained. 
Aims. To quantitatively understand the non-thermal emission process, accurate knowledge of the physical and orbital properties of 
the object is crucial. Here, we aim to investigate the orbital properties of the wide system and to constrain the inclinations of the inner 
and outer binaries, and with these the absolute masses of the system components. 

Methods. We used the PIONIER combiner at the Very Large Telescope Interferometer to obtain the very first interferometric mea- 
surements of HD 150136. We combines the interferometric observations with new and existing high-resolution spectroscopic data to 
derive the orbital solution of the outer companion in the three-dimensional space. 

Results. The wide system is clearly resolved by PIONIER, with a projected separation on the plane of the sky of about 9 milli-arcsec. 
The best-fit orbital period, eccentricity, and inclination are 8.2 yr, 0.73, and 108°. We constrain the masses of the three stars of the 
system to 63 ± 10, 40 ± 6, and 33 ± 12 Mq for the 03-3.5 V, 05.5-6 V, and 06.5-7 V components. 

Conclusions. The dynamical masses agree within errors with the evolutionary masses of the components. Future interferometric 
and spectroscopic monitoring of HD 150136 should allow one to reduce the uncertainties to a few per cent only and to accurately 
constrain the distance to the system. This makes HD 150136 an ideal system to quantitatively test evolutionary models of high-mass 
stars as well as the physics of non-thermal processes occurring in 0-type systems. 
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i1. Introduction 

HD 150136 is the brightest member of the NGC 6193 cluster in 
the Ara OB 1 association and a known non-thermal radio-emitter 
([Benaglia et al.|2006), i.e. an object where paiticles are acceler- 



ated to relativistic energies (see De Becker 2007 and references 
therein). This distinctive feature was the reason for a decade- 
long observational effort to clarify the nature and properties of 
HD 150136. Recendy, pahy et al.| ( [20T2l hereafter [Paper 1 1 
showed that HD 150136 is a triple hierarchical system consisting 
of an inner binary, with an 03V((f*))-03.5V((f^)) primary star 
and an 05.5-6V((f)) secondary star on a Pi„ x 2.67 d orbit, and a 
third physically bound 06.5-7V((f)) companion on a Pout ~ 8- to 
15-year orbit. With a total mass to be estimated around 130 Mq, 
HD 150136 is one of the most massive multiple O-star systems 
known. It is also the closest one to Earth to harbour an 03 star. 
Given the range of possible orbital periods for the outer sys- 
tem, | Mahy etaL] ( |2012[ ) estimated a probable separation on the 
plane of the sky between the inner pair and the third component 
of roughly 10 to 40 milli-arcsec (mas). Here, we report on ex- 
tended spectroscopic monitoring that allows us to obtain the first 
orbital solution of the outer system. We also report on the very 
first interferometric detection of the wide system using the Very 
Large Telescope Interferometer (VLTI). 



2. Observations and data reduction 

2.1. Spectroscopy 

We used the FEROS spectrograph mounted at the MPG/ESO 
2.2m telescope at La Silla (Chile) to obtain new high-resolution 
optical spectra of HD 150136 that supplement the FEROS data 
analysed in [Paper l] Eighteen FEROS spectra were obtained dur- 
ing a eight-night run in May 2012 (PI: Mahy), providing a con- 
tinuous wavelength coverage from 3700 to 9200A at a resolv- 
ing power of 48 000. The data were processed as described in 



Paper I In addition, we searched the ESO archives for com- 
plementary data. We retrieved six FEROS spectra from 2011 
(PI: Barba; 087.D-0946(A)), one FEROS spectrum from 2008 
(PI: Barba; 079.D-0564(B)) and one UVES spectrum from July 
2000 (PI: Roueff; 065.I-0526(A)). The latter two spectra provide 
additional constraints on the tertiary component, but not on the 
systemic velocity of the inner pair and were therefore used only 
to calculate the orbital solution of the wider pair 



* Based on observations collected at the European Southern Obser- 
vatory (ESO) under program IDs 089.D-0730 and 189.C-0644. 
** Senior Research Associate F.R.S.-FNRS 



The disentangling procedure described in Paper I for triple 
systems was applied to the 2011 and 2012 sets of spectra. 
We also reprocessed the complete data set, using our cross- 
correlation technique on all data from 1999 to 2012 to consis- 
tently measure the radial velocities (RVs) of the HD 150136 
components. The RV values coiTesponding to the new obser- 
vational epochs are given in Table 1 along with the journal of the 
201 1 and 2012 observations. 
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Fig. 1. Calibrated visibilities (left panel) and closure phases (middle panel) from PIONIER in August 2012, overlaid with the best fit binary 
model (red solid lines). The right panel shows the;^'^ map in the vicinity of the best-fit solution. 



Table 1. Journal of the new and archival spectroscopic observations. 

HJD vi V2 v^ 

-2450000 (kms-') (kms'') (kms'') 

2000 UVES observation 

1726.5548 173.5 -368.8 29.4 

2008 FEROS observation 
4658.5184 -107.5 128.2 27.0 

2011 FEROS observations 



Table 2. Interferometric best-fit measurements and Icr error-bars. 



5642.916 


-182.9 


268.4 


-6.9 


5696.790 


-208.8 


298.5 


-4.0 


5696.907 


-191.8 


269.4 


-7.7 


5697.897 


186.5 


-337.0 


-6.4 


5699.583 


-198.9 


272.9 


-6.4 


5699.588 


-191.7 


268.9 


-6.7 


2012 FEROS observations 




6048.661 


148.6 


-274.7 


-23.7 


6048.758 


120.7 


-221.2 


-27.4 


6048.899 


62.0 


-150.9 


-29.5 


6049.597 


-215.8 


306.1 


-18.0 


6049.715 


-221.1 


314.8 


-18.0 


6049.924 


-208.2 


280.1 


-22.3 


6050.666 


119.1 


-223.4 


-26.2 


6050.843 


167.7 


-310.1 


-25.7 


6051.641 


33.6 


-114.2 


-22.1 


6051.864 


-107.5 


128.2 


-21.0 


6052.646 


-192.7 


268.9 


-20.4 


6052.781 


-152.4 


219.6 


-24.9 


6052.909 


-109.9 


107.7 


-27.6 


6053.660 


188.6 


-344.1 


-27.6 


6053.778 


187.8 


-350.7 


-28.3 


6053.937 


168.4 


-308.4 


-26.1 


6054.649 


-152.4 


190.2 


-26.5 


6054.758 


-188.9 


240.8 


-23.4 



2.2. Long baseline interferometry 

Interferometric data were obtained in June and August 2012 with 
the PIONIER combiner ( |Le Bouquui et al.|2011|l and th e Aux- 
iliary Telescopes (ATs) at the VLTI ( Haguenauer et al!]|2010i l. 
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Observing date 


Parameter 


Unit 


2012-06-10 


2012-08-15 


HJD 


(-2450000) 


6088.595 


6154.573 


(h/fl+2)l.65fx 




0.24 + 0.02 


0.24 + 0.02 


6E 


(mas) 


-7.96 + 0.16 


-6.98 + 0.14 


m 


(mas) 


-4.73 + 0.09 


-5.13 + 0.10 


r 


(mas) 


9.26 + 0.19 


8.66 + 0.17 


e 


n 


239.2+1.0 


233.7 + 1.1 


^red 




1.2 


1.4 



Fringes were dispersed over three spectral channels across the 
H-band (1.50 - 1.80yum). The ATs were located in configuration 
Al-Gl-KO-Il, providing six projected baselines ranging from 30 
to 120 m and a maximum angular resolution of 2 mas in the H- 
band. Data were reduced and calibrated with the pndrs package 
( |Le Bouquin et al.]|2011| l. Calibration stars were chosen in the 
JMMC Stellar Diameters Catalog ( jLafrasseetal.|2010| l. The clo- 
sure phases and visibilities were modelled with two unresolved 
sources because the expected diameters of the individual com- 
ponents (< 0.1 mas) as well as the separation of the inner pair 
(< 0.2 mas) are largely unresolved by the longest VLTI base- 
lines. The data show no evidence that these assumpti ons may be 
wrong (xl^^ ~ 1.2). We used the LITprcFlsoftware (Tallon-Bosc 
|et al.|2 008 ) to extract the best-fit binary parameters, namely tlie 
flux ratio (assumed to be constant across the H-band) and the 
astrometric separation. The closure phases were also analysed 
independently using the method presented in |Absil et al.| ( |20lT) , 
and the results agreed excellently. 



HD 150136 was clearly resolved on both dates with a slight 
decrease in separation between June and August (Table |2]i. Fig- 
ure [T] shows the data obtained in August 2012 overlaid with the 
best-fit binary model. The final accuracy is dominated by a 2% 
uncertainty in the wavelength calibration of PIONIER. Adopt- 
ing a distance of 1.32 + 0.12 kpc ( jHerbst & HavIen|[T977| , the 
measured angular separations translate into projected distances 
of 12.2+1.1 and 11.4+1.1 AU for the June and August observa- 
tions. The Icr error-bars are dominated by the 10% uncertainty 
on the distance. 
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Table 3. Best-fit circular orbital solutions of the inner pair for the 201 1 Table 4. Best-fit LOSP RV solution of the wide system, 
and 2012 campaigns. 











Parameter 


LOSP RV solution 


Parameter 


Unit 


Primary Secondary 


Unit Inner pair Tertiary 






2011 SB2 solution 




"out 

e 


(d) 2980 + 71 

0.60 + 0.14 


Pin 


(d) 


2.67454 (fixed) 




0) 


(°) 259.2 + 5.8 


To 


(HJD- 


-2450000) 1327.136 + 0.006 




T 


(HJD-2450000) 1193 + 104 


M2/M1 




0.623 ± 0.008 




M3/M1+2 


0.31+0.07 


J 


(kms" 


1) -15.7 + 2.4 -12.1 ± 


3.2 


r 


(kms-i) -21.6 + 2.7 -16.1+4.9 


K 


(kms" 


1) 208.3 ± 2.3 334.5 + 


3.6 


K 


(kms-i) 16.7 ±4.5 53.6+14.4 


r.m.s. 


(kms" 


1) 6.2 
2012 SB2 solution 




r.m.s. 


(kms-i) 5.2 
















the RV-only 


solution, especially its eccentricity and the semi- 


/"in 


(d) 


2.67454 (fixed) 




amplitudes of the RV curves, is unfortunately limited by a lack 


To 


(HJD- 


-2450000) 1327.151+0.005 




of sampling 


around periastron. 


M2/M1 




0.630 ± 0.007 








r 


(kms" 


1) -20.4 ±2.0 -12.1 ± 


2.8 


3.2.2. Simu 


taneous RV and astrometric orbital solution 


K 


(kms 


') 214.4 + 2.5 340.4 + 


j.y 






r.m.s. 


(kms" 


1) 13.1 




The separations measured by PIONIER in June and August 2012 



3. Orbital properties 

3.1. The dose pair 

We used the Liege Orbital Solution Package^ (LOSP) to com- 
pute the orbital solutions of the inner pair during the 201 1 and 
2012 campaigns (Table [3]l. This provided a measurement of 
the systemic velocity of the inner pair at these epochs. We 
fixed the orbital period to the same value as that of [Paper l) i.e. 
2.67454 days. For both campaigns, the mass ratio, the semi- 
amplitudes, and the projected semi-major axis are similar to 
those found in Paper I. The difference between the systemic ve- 
locities of the piimary and of the secondary components is larger 
for the 2012 solution but, globally, both orbital solutions of Ta- 
ble[3]ag ree with each other and with the solution determined in 
[Paper I| 

3.2. The third companion 

3.2.1. Spectroscopic orbital solution 

We determined the systemic velocities of the inner system during 
the individual short-term campaigns with more than one spec- 
trum since 1999. These values are put in perspective with the 
RVs of the tertiary component in Fig. l2] There is a clear, anti- 
correlated periodic motion of the systemic velocity of the close 
binary system and of the tertiary component. We performed a 
Fourier analysis on the tertiary RVs using the Heck-Manfroid- 
Mersch method Peck etal.|l 985 1 refined by [Gossetetal.|2001 ). 
The highest peak in the periodogram indicates a period of about 
3335 ± 260 days. We adopted this value as our initial guess and 
used LOSP to refine the period estimate. We computed an SBl 
orbital solution from the RVs of the tertiary alone and an SB2 
solution using the tertiary RVs and the systemic velocities of the 
inner pair ijn). The two solutions are in excellent agreement 
and the SB2 RV solution is given in Table |4] The reliability of 



' www.jrmnc.fr/litpro 

^ LOSP is developed and maintained by H. Sana and is available at 
http://www.science.uva.nl/~hsana/losp.html The algorithm is based on 
the generalization of the SB 1 met hod of, Wolfe et al.H1967| l to the SB2 
case along the lines described in |Rauw et al.| ( |2000^ and [Sana et al.| 
1 2006} . 



are relatively small and indicate an eccentricity at the upper end 
of the confidence interval of the RV-only solution of Table |4] 
We therefore proceeded to simultaneously adjust the RV and as- 
trometiic measurements of the outer pair. We minimized the 
squared differences between the measurements and the model 



X 



0-12 j ^\ a 
y;-^ ( 6E - 6E"'°'^\^ 



y(n2^y]^\ 



..mod y^ 



OSe 






(1) 



using a Levenberg-Marquardt method, adopting the RV-only so- 
lution of Table |4] as a starting point. Based on the residuals of 
the tertiary and the systemic RVs around the best-fit RV curves, 
we adopted 0-3 = 3.4 and 0-12 - 6.8 kms"' as typical uncertain- 
ties on ut, and y^. The three-dimensional orbital solution con- 
verges towards a higher eccentricity (although still within error- 
bars), hence towards larger RV curve semi-amplitudes, than the 
RV-only solution. All other parameters remain unchanged. The 
best-fit combined solution is given in Table l5] and is shown in 
Fig.|2] 

To estimate the uncertainties, we performed Monte-Carlo 
(MC) simulations. We randomly drew input RVs and astromet- 
ric positions around the best-fit orbital solution at epochs coiTe- 
sponding to our observing dates. We used normal distributions 
with standard deviations corresponding to the Icr measurement 
uncertainties. We also included the uncertainty on the distance 
by drawing the distance from a normal distribution centered on 
1320 pc and with a standard deviation of 120 pc. We recomputed 
the best-fit orbital solution and inclination using 1000 simulated 
data sets and we constructed the distributions of the output pa- 
rameters. The medians of the distributions match the best-fit val- 
ues of Table |5] very well. We used the distances between the me- 
dian and the 0.16 and 0.84 percentiles as uncertainty estimates. 
The upper and lower error-bars agree within 10% for all param- 
eters; Table |5]provides the average of both values. 

4. Discussion 

4. 1 . Physicai properties 

Using the best-fit three-dimensional orbit, absolute mass esti- 
mates are M3 = 33 ± 12 M^ and M1+2 = 102 ± 16 M^ for 
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Table 5. Best-fit simultaneous RV and astrometric orbital solution of 
the wide system. The corresponding RV curves and projected orbit on 
the plane of sky are displayed in Fig. [2] 



Parameter Unit 



Combined solution 
Inner pair Tertiary 



"out 


(d) 


e 




(0 


(°) 


T 


(HJD-2450000) 


M3/M,+2 




'out 


(°) 


Q 


(°) 


J 


(kms-i) 


K 


(kms-i) 


a 


(A.U.) 


M 


(Mo) 


-^red 





3008 + 54 
0.73 ± 0.05 

250.7 + 2.9 
1241 + 38 
0.32 ± 0.08 

107.9 ± 3.0 
-60 ± 10 



-19.9 ±2.2 

25.8 ±7.6 

5.1 ±1.2 

102 ±16 



-21.5 ±2.0 
79.8 ±11 
15.8 ±0.6 
33 ±12 



0.94 



the third companion and the total mass of the inner pair, respec- 
tively. Combining the total absolute mass of the inner pair with 
the spectroscopic minimum masses of [Paper I] (table 4), we con- 
strain the inclination of the inner binary to values of 49.6 ±3.6°. 
This value is compatible with the absence of eclipse in the in- 
ner system. Best masses for the primary and secondary are thus 
Ml = 62.6 ± 10.0 and Mj = 39.5 ± 6.3 Mq. The agreement 
with the expected masses given the estimated spectral types is 
remarkably good ( Martins et al.|20()5] l. The dynamical masses 
also agree within errors with the evolutionary masses obtained 
in |Paper Ij 70.6+^\•^ 36.2+f° and 27.0+i^° M^ for components 1, 
2, and i, respectively. The measured flux-ratio in the H-band, 
(/3//i+2)i.65/i = 0.24, also agrees well with the expected H-band 
flux-ratio of 0.26 given the component spectral types (Martins 
|& Plez||2006[ ), which confirms the main-sequence nature of the 
tertiary object (see discussion in Paper I| i. 



4.2. Non-thermal emission 

The present study allows us to clarify to some extent the consid- 
erations on the origin of the non-thermal radio emission. As ar- 
gued in Paper I, the synchrotron radio emission probably comes 
from the colliding winds in the wide orbit. The radius of the 
photosphere (t - I) for the primary (whose wind dominates the 
free-free opacity in the system) is expected to be at most 850 Rq 
(at A - 20 cm), strongly suggesting that the stagnation point of 
the collision zone must be located farther away (see Paper I). For 
the tertiary star, we assume two typical va lues for the mass loss 



rate: M^; = 10" Moyr~ (classical value, Muijres et al. 2012 i 
and Mu,„ - 10"^ Mq yr"' (more representative of the weak-wind 
case). We also adopt a terminal velocity of 2500 kms"', corre 
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Fig. 2. Main panel: Evolution of the third component RVs (triangles) 
and of the systemic velocity of the O+O inner system (squares) as a 
function of time. The best-fit RV curves are overlaid (Table [5|. Insert: 
PIONIER astrometric points and projection of the best-fit relative orbit 
on the plane of the sky. The dashed line shows the line of nodes. 



the assumption about the nature of the stellar wind (classical or 
weak). 

In the classical- wind case, however, the stagnation point very 
close to periastron could be located below the radio photosphere 
at longer wavelengths (13 and 20 cm) where HD 150136 is sig- 
nificantly detected, but the observations reported by [Benaglia] 



et al. (20061 were performed in December 2003, which is far 
from periastron passage according to our ephemeris. In the 
weak-wind case, the stagnation point is always located signifi- 
cantly outside the photosphere of the dense primary wind, even 
at periastron, and at all wavelengths. These facts are in agree- 
ment with the rather high flux densities measured between 3 and 
20 cm, as reported by "Benaglia et al."('2006). 

As discussed in De Becker et al. ( 2012) for a similar system, 
the variability in the synchrotron radio flux from a colliding- 
wind binary most likely comes from two changing factors: stel- 
lar separation and free-free absorption. Following the discussion 
above, absorption is unlikely to dominate. Any temporal change 
of the HD 150136 flux densities would therefore be mainly at- 
tributed to the varying stellar separation along the eccentric wide 
orbit. One may therefore expect the radio flux density to reach its 
maximum very close to periastron, with a minimum at apastron. 

A radio monitoring of the wide orbit in HD 150136 (so far 
non-existent) is required to validate this scenario and to achieve a 
more detailed description of the relevant non-thermal processes. 
Given the expected location of the stagnation point, the wind 
interaction zone may be resolvable from the stars by very large 
baseline radio interferometric facilities. Imaging the wind-wind 
collision at several epochs would help to distinguish between a 
weak-wind, as possibly suggested by the CMFGEN analysis in 



sponding to 2.6 x ^esc and Wesc 
( [Muijres etal.|2012| l. 

On these bases, we estimate that, in the classical-wind case, 
the stagnation point is located at about 260 and 16307?© from 
the tertiary at periastron and at apastron. In the weak-wind case, 
these distances reduce to about 30 and 200 Rq. Using the same 
approach as described in Paper I, we find that the extension of 
the radio photosphere of the tertiary star is shorter than the dis- 
tance between the tertiary component and the stagnation point 
of the colliding winds, whatever the wavelength and whatever 
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970 km s ' for 06 5-7 V stars IP^P^^ ^1 ^"^ ^ normal wind for the tertiary star. 



5. Summary 

We reported the very first interferometric detection of the outer 
companion in the hierarchical triple system HD 150136. Com- 
bining the interferometric measurements, new and archival spec- 
troscopy with data from Paper I we obtained the first three- 



dimensional orbital solution of the wider system. The best-fit so- 
lution indicates a 8.2-yr period and a high eccentricity (e > 0.7). 
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The accuracy of the PIONIER interferometric measurements al- 
lowed us to constrain the inclination of the outer orbit, and from 
there on, of the inner pair to within a few degrees only. 

We constrained the masses of the three stars of the system 
to 63, 40, and 33 Mq for the 03-3.5 V, 05.5-6 V, and 06.5-7 V 
components. In particular, this is the first direct measurement of 
the mass of an early main-sequence O star We showed that the 
obtained dynamical masses agree within errors with the evolu- 



tionary masses estimated in Paper I 



Although relative error-bars remain at the 15%-level for 
masses of the inner pair components and at the 30%-level for 
the tertiary mass, spectroscopic observations around the next pe- 
riastron passage in 2015 and further interferometric monitoring 
over a significant fraction of the orbit may provide direct mass 
measurements with uncertainties of a couple of per cent only. It 
will also allow for an accurate independent measurement of the 
distance to the system. This constitutes a must-do to accurately 
test high-mass star evolutionary models and to investigate the 
reality and origin of the mass discrepancy in more details (e.g. 
IHerrero et al.fT992l|Weidner & Vink|2010l l. 

As for other non-thermal emitters (HD 93250, [Sana et aL] 
[20TT] HD 167971, De Bec ker erai:][2nT2| , the present study 
also shows that long baseline interferometry is ideally suited to 
resolve the components of O-type non- thermal radio emitters. 
The significant advances on the three-dimensional description of 
the orbits in HD 150136 is a solid base for future investigations 
to understand the non-thermal emission and particle acceleration 
processes at work in massive multiple systems. 
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